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THE EFFECT OF FINITE THRUST AND HEATING CONSTRAINTS ON THE 
SYNERGETIC PLANE CHANGE MANEUVER FOR A SPACE SHUTTLE 
ORBITER-CLASS VEHICLE 


INTRODUCTION 


In the synergetic plane change maneuver, propulsive and aerodynamic forces are 
used in conjunction to achieve a rotation of the orbital plane by dipping into a planetary 
atmosphere using the lift to effect part of the plane change, the thrust providing the rest and 
adding the lost energy due to drag. Several special maneuvers have been discussed in t'he 
literature since 1962 [1-18]. A survey over these has been given in Reference 19. In all 
investigations, the total maneuver has been broken down into deboost, descent, 
atmospheric/ thrusting flight, ascent, and injection into the new orbit. To narrow the field of 
parameters, mainly plane changes of circular orbits have been considered. The phases 
mentioned above were investigated separately and later pieced together, yielding only rough 
estimates of the total performance achievable. 

In References 19 and 20, a variational formulation for the atmospheric skip phase 
was given and analytical approximate and numerical results for unconstrained maneuvers 
and the effect of a normal acceleration constraint were presented. The deboost phase and 
the ascent bum phase were considered to be impulsive which together with information 
from the aerodynamic glide phase allowed an approximate determination of favorable 
interface parameter between these phases. Especially, ascent thmst angles and thrusting 
altitudes were obtained from simple function minimizations. Kinetic heating considerations 
had not been considered since the estimation of a lower limit for the AV requirement was 
the primary purpose of the investigation. The resulting aerodynamic skips showed high 
heating and acceleration peaks, especially for larger plane change angles which desired 
steep entry flight path angles. 

For a realistic maneuver with state-of-the-art materials, a heating constraint has to be 
considered which is more severe than the normal acceleration constraints considered in 
References 19 and 20. This type of performance reduction caused by the vehicle skin 
temperature limitation has been investigated by Clauss and Yeatman [12] for a minor circle 
glide turn and a steady state aerocruise maneuver. They concluded that the temperature 
constraint is more detrimental to the glide than to the aerocruise maneuver. 

The objective of this report is to investigate the effect of a heating limit on a finite 
thrust optimal plane change maneuver for a Space Shuttle orbiter-type vehicle. Only the 
central portion of the maneuver (Fig. 1) which lies in the dense atmosphere is considered; 
this, however, is a full variational formulation. The remaining mass fraction after the burn 
phases is being maximized. Numerical result for a limited range of parameters will be given. 
They were obtained partly by a refined gradient optimization program and partly by a 
multiple shooting algorithm to solve the boundary value problem resulting from the 
application of the maximum principle. 



MATHEMATICAL MODEL 


To keep the computational work load in this exploratory study low, a relatively 
simple process model has been adopted which is described in this section. 


Planet Characteristics 

The planet is assumed to be a nonrotating, sphere with a stationary atmosphere 
independent of geographical coordinates. No winds are assumed and the density altitude 
relationship is exponential 


p = p Q exp (-0h) (1) 

# 

where 0 is the inverse scale height 0 = 1. 0/6.9 [km" 1 ] and p Q is the density at sea level 
p Q = 1.54 [kg/m 3 ] for earth. 

An inverse square gravity law 


g=r/r 2 

(2) 

r=R 0 + h 

(3) 


with sea level radius R Q = 6371.2 km for earth and T = gravitational constant times 
planet mass (3.986 X 1 0 14 m 3 /s 2 for earth) were considered a good approximation. 

The neglection of the planet rotation results in large errors in the heating equations 
for equatorial entry, where, for earth the velocity increment caused by rotation is 6 
percent of the satellite velocity resulting in about 20-percent heating rate error. Therefore, 
the results presented correspond more to a meridional entry in higher latitudes. Viscous 
interaction effects at high altitudes are neglected, leading to a constant hypersonic drag 
polar for the vehicle. 


Vehicle Model 

The drag polar was chosen as a constant medium one for the altitude range of the 
maneuver and represented by 


C D = C Do + k C L " 


(4) 
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where Cjj 0 is the zero lift drag coefficient and n is the power of the lift-dependent drag 
(actual values: = 0.04, k = 1.0, and n = 1.86, resulting in a maximum L/D of 2.22 for 

^LE = 01 92). The considerable influence of viscous interaction effects on the optimal 

aerodynamic control for the unconstrained synergetic maneuver has been shown in 
References 19 and 20. The effect on the performance, however, was only minor. 

The second important aerodynamic parameter is the wing loading. Low values lead 
to higher skipping altitudes and lower heating loads while the vehicle has to be larger for the 
same payload. A value of m 0 /F = 275 [kg/m 2 ] was chosen as nominal; the effect of 

increasing it to 300 kg/m 2 was determined by one sensitivity test run. No aerodynamic side 
forces were considered assuming that the sideslip angle is always zero. 

A chemical propulsion system with unconstrained thrust direction was assumed 
without atmospheric back-pressure losses 


T 

nu 


= b - LL 


(5) 


The exit velocity U e was chosen to be U e = 4.36 [km/s] and the relative mass flow rate 
b ^ 0.002 [s' 1 ] , resulting in an initial thrust acceleration of about 1 g Q . 


Heating Constraint 

The heating constraint for a reentry vehicle is one of the driving factors in trajectory 
shaping. Medium- (like the proposed Space Shuttle orbiter) and high-L/D vehicles have a 
thermal protection system (TPS) consisting of reradiative and ablative elements at different 
parts of the vehicle [21 ] . The higher the L/D, the more reradiative the TPS is going to be. If 
the time constants for the outer skin to heat up to the design limit is small compared to the 
rate at which the vehicle state (especially altitude, angle of attack, and velocity) changes, 
then a quasisteady approximation is valid. In this case, the heating constraint can be 
formulated as an algebraic equation 0(v,h,a) = 0, where v is the aerodynamic 
velocity,, h is the altitude (representing air density), and a is the angle of attack of the 
vehicle. It is assumed that the sideslip angle 0 is kept small so that its influence can be 
neglected. 

The data underlying the present representation are taken from Reference 22 for a 
limit temperature of 1093°C (2000° F). Figure 2a shows a qualitative picture of typical 
altitude constraints caused by dynamic pressure and kinetic heating. The effects of both 
velocity and lift coefficient are seen to be appreciable for the heating constraint. The 
heating constraint is more severe than the dynamic pressure boundary for hypersonic speeds 
(v > 2.5 to 4.5 km/sec, depending on angle of attack). 

In the trajectory optimization program, the control Cjj_j on a boundary arc is 
required as a function of the state variables altitude and velocity. 
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The data are given in the form [22] 


0(h, v, a) = 0 (6) 

Since the use of the lift coefficient instead of angle of attack a reduces the 

computational workload to represent the aerodynamic characteristics of the vehicle, the 
relation CjTa , h , v) is used to eliminate a . The transformed equation (6) is then 

represented in the form 


C LH -B i H i + AC LH 


i = 1, 5 


(7) 


where 


B i = G iJ ^ (j - 1) , j = 1, 4 
H, = b 2 h 2 /v 


H 2 = bh/v - H! 


( 8 ) 


H 3 = 1 - bh/v - H 2 


H 4 = v/(bh) - 2 + bh/v - H 3 

H s = v 2 /(bh) 2 - 3v/(bh) + 3 - bh/v - H 4 


Clh is an adjustment parameter. 

Gj j is a coefficient matrix as 

.110717 

.834519 , 

1.213679 

, -1.060833 

-.672677 

2.73417 

-.864369 

, -12.100000 

.812241 

2.337815 , 

10.316280 

, 22.974860 

-3.151267 

-13.621310 , 

-40.485500 

, -57.833330 

2.368095 

19.073400 , 

69.86905 

, 127.777778 


( 9 ) 
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and h is a coordinate transformed altitude 


n = h/50 [km] - 1 


( 10 ) 


to keep the coefficients closer to 1. The constant b = 0.095 serves the same purpose for the 
velocity. This form has been obtained by intuition and trial. It is a bivariate polynomial of 
third order in the altitude and fourth order in velocity. With 

C r t G-DGi/hCi- 2 ) , (11) 

j=l 

the partial derivatives may be written 


3C L h/ 3v = Bja/av (Hj) 

( 12 ) 

3C Lh /3h = Cj Hj + Bj 3/3h (Hj) 

(13) 


From the form ( 6 ), one obtains the partial for constant 


3h/3v 


C LH 



(14) 


The function generated as above is shown in Figure 2b. 

The accuracy of the approximation is indicated in Figure 2b by dots representing 
the original input data into the curve-fit procedure. These points had been obtained from 
Reference 22 through a crossplot. Except for the 70-km curve, the approximation is good 
considering the uncertainty in the original data. For maximum range and minimum 
energy-loss trajectories, which result in lift coefficients in the vicinity of maximum L/D 
(around C ^ = 0.2), the 70-km curve will not be needed. Aside from that, the deviation is 

on the safe side. Below v = 2.5 [km/s] , the constraint equation is not being used. 


Equations of Motion 

For the coordinate system, a flight-path-oriented axis system is chosen (Fig. 
3). The x-axis is tangential to the trajectory, positive in the flight direction. It is inclined 
relative to the horizon by the flight-path angle 7 , positive upward. The heading 
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angle x designates the orientation of the horizontal velocity component relative to the 
initial time t Q , positive from north over east. The range angle in the original orbit plane 

is < 3 > and A (lateral range) is normal to it. 

Five controls exist: two associated with the aerodynamic forces and three with the 
thrust forces. The magnitude of the lift coefficient and the lift bank angle p > positive 

from the vertical for increasing heading, determine directly the plane changing force 
component sin ju Fq . The lift coefficient also controls the drag according to equation 

( 4 ). The three thrust controls are the thrust magnitude, for fixed jet velocity [equation ( 5 )] 
governed by the mass flow rate b = rh/m 0 , the plane-changing thrust angle e , measured 

from the vertical plane containing the velocity vector, and the thrust angle a in the vertical 
plane, measured from the velocity vector. This leads to the following system of equations of 
motion: 


Inertial 


Atmospheric Propulsive 


v = - g sin 7 


-D/m 


— cos e cos a 
m 


v . . , L sin ju 

— cos 7 cos x tan A + 

r mv cos 7 


T sin e 


mv cos 7 


(?-*) 


cos 7 


mv 


sin n + 


mv 


cos e sin o 


<S> = - 


v cos 7 


r cos A 


cos x 


05 ) 


A 


— cos 7 sin x 
r 


Kinematic relations 


h = v sin 7 


m = - b 


where 
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L 

D 




are the aerodynamic force components, lift and drag. 


Boundary Conditions 

A complete synergetic plane change maneuver only has the specifications of the 
initial and final orbit as boundary conditions. By removing the legs of the total maneuver 
that are outside of the denser atmosphere, as is done here to reduce computational 
difficulties, somewhat arbitrary boundary conditions are introduced. As will be shown in the 
discussion of the results, some complications may be encountered because of this procedure. 

Using engineering judgement, the range of entry and exit parameters can be 
narrowed considerably. The “edge” of the denser atmosphere may be defined as the altitude 
in which the total aerodynamic force grows to a threshold percentage of the gravitational 
force. 


In this study the altitude, h = 80 km , has been selected even though at this level 
heating is already significant. Velocity losses caused by drag above this level can be 
estimated analytically [ 19] . 

From the qualitative picture of the heating constraint in Figure 2a, it may be seen 
that minimum altitudes during the skip cannot be expected below 30 km. The entry 
velocity for a return from a 550 km (300 nm) circular earth orbit is around 7.95 km/s. 
From these numbers, an estimation of favorable entry flight path angles may be obtained. 
From energy and angular momentum relations for elliptic vacuum orbits using small angle 
approximations (sin 7 e « 7 e , cos 7 e * 1 - 7 e 2 /2) , one obtains 



(16) 


where v $e is the satellite velocity in the entry altitude and Tp is the perigee radius. The 
entry angles for v e = 7.95 km/s and entry altitude of 80 and 95 km corresponding to 
perigee altitudes of 30 to 60 km are given in Table 1. 

It is seen that for this entry velocity, the flight path angle at 80 km is about 0.25 deg 
shallower than at 95-km altitude. To just hit the heating constraint with the vacuum perigee 
altitude at maximum L/D, the entry flight path angle should be around -0.7 to -0.8 deg. This 
small flight path angle, however, results in a long descent arc, long maneuver time and larger 
total heat inputs. Therefore, a nominal entry angle of -1.25 deg at 80 km has been chosen. 
This results in conservative trajectories since both the descent and the ascent propellant 
requirements are somewhat higher than that required for the maneuver. The effect of entry 
flight path angle on the atmospheric flight arc will be discussed later. 


7 



TABLE 1. ENTRY FLIGHT PATH ANGLES IN DEGREES FOR v e = 7.95 km/s 


^\hp/km 





h e /km 

30 

40 

50 

60 

95 

-1.393 

-1.266 

-1.13 

-0.983 

80 

-1.151 

-1.015 

-0.866 

-0.697 


The exit boundary conditions for the heading angle x and the lateral 
range A determine the plane change angle i . Combining the inertial term of the second 
equation with the fifth equation in equation ( 1 5) yields 


cos x * cos A * = cos x cos A 


(17) 


after integration. For A* = 0 , the heading angle x* is identical to the plane change 
angle i . It is assumed that on the upgoing leg beyond the specified boundary values of Vf = 

7.95 km/s and hj- = 80 km no further plane change takes place. The relation between 
plane change angle i , heading angle at exit Xf » and the lateral range Af in the parameter 
range of interest here is plotted in Figure 4. 

In the numerical iterations,' the final heading angle Xf had been specified, and the 

lateral range was left free to assume the most favorable value. The full dots in Figure 4 show 
that the plane change achieved is approximately a fixed percentage (« 1 .4 percent) higher 
than the specified heading for the unconstrained trajectories. In the constrained cases (light 
circles), the percentage increases with specified final heading angle. When comparing both 
cases, this difference has to be considered. 

Resulting from similar considerations with respect to maneuver range and total heat 
transfer as for the entry, the exit flight path angle has been selected to be 7 f = 1.5 and 2 
deg, respectively (see also Fig. 12 of Ref. 20). The specification of 7 f may have a strong 

influence on the optimal bum time as will be discussed later. Small exit angles may result in 
a split of the ascent burn phase. 

The approach taken results in the following set of boundary conditions for the 
system of differential equations (15): 
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Variable 

Initial Conditions 

Final Conditions 

v/km/s 

7.95 

7.95 

x/° 

0 

Specified 

y/° 

Specified 

Specified 

A/° 

0 

Free 

h/km 

80(95) 

80(95) 

m/m 0 

1 

To be maximal 


The downrange angle © does not enter the differential equations and was therefore 
omitted. 


REDUCTION TO A BOUNDARY VALUE PROBLEM BY 
APPLICATION OF THE MAXIMUM PRINCIPLE 


The nonlinear optimal control problem with six state and five control variables is 
reduced to a boundary value problem, which has to be solved numerically, by the applica- 
tion of the maximum principle. 


The Hamiltonian function is from equations (1), (5), and (15): 

H - y cos t ^ X T (l - (X<2> -* x sin x ) + X A sin x + sin 7 (vX h - gX 


p Po 


2m 


+ ve“0 h | Cr ( X 

J-A X cos 7 


+Xt cosl j " vX » c d 


+ b] 

K 

(vX., cos a + X., sin a ) 

x x • 1 
cos e + — sin e 

-x m ) 

1 

| mv 

V v 7 / 

cos 7 



(18) 


The auxiliary (costate) variables X are given by the differential equations 


X = -9H/3X 


( 19 ) 


which are derived in the appendix. 
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Determination of the Optimal Controls 

Since the payoff quantity final mass is to be maximized, the Hamiltonian has to be 
minimal with respect to all controls 


U T = (C L , n , b , a , e) 


Extremal Bank Angle. The expression ^X^/cos y ^ sin p + X^ cos ju has its 
minimum -w„ with respect to u 



( 20 ) 


for 


sin q* = X x /(-w a cos 7) , cos n* = -\ y /w a ■ (21) 

Extremal Lift Coefficient. The square-bracket expression in the second line 6f 
equation (18) contains all lift coefficient terms. Using equation (4) for the drag polar 
yields 3H/3Cy = 0 for 


C L *= [-w a /(vX v kn)] 1//(n_1) 


( 22 ) 


Extremal Thrust Magnitude. In the third line of equation ( 1 8), the mass flow b 
appears linearly so that the expression in the outer bracket acts as a switching function 


Sy = U e (-wj)/(mv) - X 


m 


(23) 


Depending on the size of Sy , the following values for control b have to be chosen 
(for wy see below) : 


Sy X 0 
Sy = 0 
Sy 0 


b = 0 

Intermediate thrust arc (singular case) 
b - ^max 


(24) 
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Extremal Thrust Direction in Vertical Plane. The innermost bracket in the third line 
of equation (18) has its minimum with respect to a for 

sin a* = -A /w, ; cos a* = -v\ v /wj (25) 

with 

-w, =-[(vX v ) 2 +V ] 1/2 • (26) 

Extremal Plane Changing Thrust Angle. With equation (26), the middle bracket in 
the last row of equation (18) may be written as 

-Wi cos e + (x^/cos y') sine , 

which has its minimum with respect to e for 

sin e* = A x /(-wj cos 7 ) ; cose* = W!/wj (27) 

with 

-wj = - jw , 2 + (A x /cos t) 2 J ^ (28) 

From equations (25) and (27) follows 
cos e* cos a* = -v\ v /wj 


cos e* sin a* = -X^/wj . (29) 

Now, the extremal controls are given at each instant as a function of the state and 
costate variables. 
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For each unspecified boundary value of the state variables there is a relation given 
from the transversality conditions for the corresponding boundary values of the costate 
variables so that in total 2n boundary values specify the solution to the 2n-system of 
differential equations (15) and (19). The requirement Hf = 0 determines the final time. 


Constraint Arcs 

On constraint arcs involving one control directly, the constraint equation 


C(X, u) = 0 


(30) 


determines the control as a function of the state [23]. For the constraint to remain 
satisfied, the following relation has to hold 


3C . . 9C . n 

-r — 5 x + — 5 u = 0 

3x 3u 


(31) 


This yields 




ac 

ax 


sx 


(32) 


which when introduced into the linearized perturbation equation 


5X= 6X+ -P- 5u 

ax 3u 


(33) 


leads to 


5X= 


_3f 

ax 


l/ac) 

u \3u ) 


ac 

ax 


6 X 


(34) 


The costate (adjoint) differential equations therefore have additional terms on constraint 
arcs, where the control is being determined from the constraint equation 


X = - 


af 

ax 


_ II /ML) 1 

au \au/ 


ac 

ax 


(35) 
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For the heating constraint (6) which may be written in the form 


C - C L - Cjjj(h, v) < 0 


( 36 ) 


the following terms are obtained: 


3C = 
3X 

/T_ 


ac 


LH 


3v 


0 0 0 


ac 


LH 


ah 


where X 1 = (v, x, 7, A, h, m), 


(37) 


ac 

3Ct 


= l 


(38) 


and from equation (15) [X = f(X,u)] results 


II 

3u 


ac 

3u 


-l 


ac 

ax 


= _ Ip 


2m 


-v 


aC D aC LH 

ac. 


3v 


sin n aC LH 


cos 7 3v 


ac 


COS n 


LH 


3v 


0 

0 

0 


0 0 0 -V 


0 0 0 


aC D aC LH 


ac, ah 


sin ju aC LH 


cost 9h 


0 0 0 cos n 

0 0 0 

0 0 0 

0 0 0 


ac 


LH 


ah 


o 

o 

o 


(39) 


This yields, as additional terms in the costate equations on constraint arcs, 


\ = * v + 
c v uc 


ac 


LH K 


3v 


(40) 
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(41) 


X 


he 



ac 


LH 


ah 


K 


where 



ac 


w a + v V 3C 


D\ Fp 


2m 


SOLUTION OF THE BOUNDARY VALUE PROBLEM 
BY MULTIPLE SHOOTING 


(42) 


The majority of the numerical solutions discussed in the next section were generated 
by a multiple shooting algorithm [24], Initial gliding solutions were obtained by a min-H 
gradient program. Then, by using continuation methods in several parameters, the 
glide-thrust trajectories were generated with a general multiple shooting subroutine due to 
Bulirsch [25]. This isolation scheme has convergence properties ranging from those of 
simple shooting to those of the generalized Newton-Raphson algorithm depending on the 
number of segments used [25, 26] . It may be viewed as a compromise for combining low 
storage requirements and easy adaptation to new problems of the simple shooting with the 
possibility of generating iterates not from just one set of variables at one point in time. By 
the introduction of intermediate grid points, the sensitivity of variables at one boundary 
with respect to slight changes of variables at the other boundary may be reduced 
significantly. Numerical errors propagate only over each segment. Figure 5 shows 
qualitatively the gain to be expected. Only one variable is displayed. The upper two curves 
for simple shooting show that the computation of the Jacobian matrix by numerical 
differentiation is not achieved since the trajectory based on poor guesses, for the missing 
initial conditions blows up, a case easily encountered in reentry. With multiple shooting, the 
trajectory may be restarted at intermediate points using new guesses. Jacobian matrices are 
formed for each segment. The resulting iteration scheme, based on reducing all 
discontinuities at internal grid points to zero, leads to a system of linear algebraic equations 
which can always be reduced to the dimension of the simple shooting problem. The total 
Jacobian matrix linking the changes of the boundary values at each end to each other is 
formed as the product of all segmental Jacobians. Through the intermediate grid points, the 
path along which the total Jacobian is formed can be much better controlled. The price to 
pay, however, is that for interior segments the Jacobian matrices required are full square 
matrices, whereas in simple shooting and, of course, for the first segment in multiple 
shooting, only rectangular matrices diminished by the known initial conditions are required. 

Three to five segments, depending on the length of the maneuver, were selected. In 
regions where aerodynamic forces are significant, segment lengths of 100 to 150 s seem 
appropriate. On the descending or ascending leg much larger segments could be handled. 
Figure 6 shows a typical example of the segmentation and iteration. The initial guesses for 
the second and third segments were relatively poor resulting in large jumps at the 
intermediate grid points (Figs. 6a and 6b). They are reduced considerably during the first 
iteration. After convergence, the curves are smooth and the altitude-velocity phase plot ends 
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at the same point where it started (h = 95 km, v g = 7.95 km/s, not shown; for similar plots 

see Fig. 9b). The time histories (Fig. 6c) show that the intermediate grid points are at the 
beginning, the middle, and the end of the high dynamic pressure region of the trajectory. 
For vacuum trajectories simple shooting often is a good iteration scheme as can be inferred 
from the long, low density descending and ascending portions. Without the intermediate 
grid points, however, convergence will often not be achieved. In the iteration scheme used, 
the time was normalized to the range 0 4 1 . Maneuver time adjustment was derived from the 
requirement that the Hamiltonian function at the final time has to be zero for open final 
time. This adjustment shows up as a shift of the grid points in an absolute time scale (Fig. 
6c). 


DISCUSSION OF RESULTS 


Unconstrained Finite Thrust Solutions 

Figure 7 shows the altitude time histories and the optimal controls for final heading 
angles up to 40 deg. In Figure 9, phase plots in comparison to constrained trajectories are 
given. Figure 7a indicates the increasing penetration depth into the atmosphere with 
increasing plane change. The main thrusting arcs are on the upgoing legs of the skip, and 
center around an altitude of 60 to 65 km. For the larger plane change angles, the thrusting 
arcs extend to the right boundary. For certain parameter combinations, the iteration calls 
for a coast arc between two thrust arcs on the upgoing leg as indicated by the switching 
function for Xf = 35 deg in Figure 7b around 500 s. For slightly smaller heading changes, 

this coast arc would show up. As will be seen later, this upward bending of the switching 
function seems to be induced by the specification of too small a final flight path angle. 
From the thrust angle a in the vertical plane (Fig. 7c), it may be concluded that the 
upward accelerating component of the first thrust arc helps to reduce drag losses by 
steepening the flight path angle while the second (short one) at the end with a downward 
accelerating component (Fig. 7c) helps satisfy the specified end conditions. A similar 
behavior occurs for Xf = 12.5 deg and 7f = 1.5 deg. In these cases, the convergence 

behavior is very poor. Increasing 7 f to 2 deg at Xf = 12.5 deg moves the occurence of the 
second burn arc to higher final heading angles. 

From Figure 7b, it is seen that for small plane change angles (xf X7.5 deg) an initial 
burn phase appears for the entry angle 7 e = -1.25 deg chosen. It both decreases the 

flight path angle and reduces the wingloading. The effect on the flight path angle is 
two-fold: a direct influence caused by the upward thrust component (see circles at t = 0 in 
Fig. 7c) and an increased centrifugal force caused by the velocity increase. For smaller entry 
flight path angles, this thrust phase vanishes. With increasing plane change, which leads to 
increasingly steeper optimal entry flight path angles [20, Fig. 7] in the unconstrained 
case, the initial value of the switching function for the fixed entry flight path angle 
considered here increases, indicating that initial bum arcs become less favorable (Fig. 7b). 
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The plane changing thrust angle e in the upper part of Figure 7c varies around 
23 deg with its highest value at the beginning of the thrust phase when the velocity 
vector to be turned is smallest and decreasing with increasing time and velocity. This 
optimal thrust angle is well approximated by the analytical solution of References 19 and 
20 or the even simpler result e*= tan" 1 (D/L) which yields for L/D = 2.22 used here a 
value of e* = 24.25 deg. 

The optimal aerodynamic controls are shown in Figures 7d and 7e. The bank 
angle time histories are similar to those for gliding flight with maximum exit velocity at a 
given heading and flight path angle [19,20], The lift coefficient time histories have 
characteristic kinks at thrust onset and cutoff. During thrust periods they decrease 
monotonically. The lift coefficient is always in the vicinity of the value for maximum 
L/D. 


Synergetic Plane Change with Heating Constraint 

Figure 8 shows results for the nominal parameter set chosen : entry and exit 
altitude = 80 km, velocity = 7.95 km/s (v se = 1 .01 ), entry flight path angle 7^ = -1 .25 deg, 

initial wing loading m Q /F = 275 (kg/m 2 ), initial thrust acceleration about 1 g Q . 

Comparing Figure 8a to Figure 7a, it is seen that the minimum altitude decreases much 
slower with increasing plane change because of the heating constraint. For Xj- = 17 deg, 

the minimum altitude is 60.45 km compared to 52.7 km in the unconstrained case. The 
maximum dynamic pressure is reduced to about 6700 N/m 2 (one third of the 
unconstrained maximum) and the maneuver time is increased by about 65 percent. The 
maximum normal acceleration is approximately 0.43 g Q compared to about 1 .6 g Q in the 

case with no heating constraint. 

A short initial burn arc appears up to higher plane change angles than in the 
unconstrained cases; for the entry parameters chosen, it vanishes for Xj-^13deg (Fig. 

8b), where the thrust switching function becomes negative initially. If a steeper entry 
flight path angle would be chosen, the initial burn phase would increase, in time and show 
up for larger plane change angles. The general shape of the thrust switching function is as 
follows: It starts with a negative slope, passes through a minimum and a maximum and 
may have a second minimum (not shown), depending on the specified final flight path 
angle 7^.. This may occur when 7^. is smaller than the selected value of 2 deg. From 

Figure 8b it is seen that the slope of the switching function at zero crossing between the 
first minimum and the maximum, i.e., at the beginning of the ascending leg of the 
trajectory, becomes smaller for increasing plane change angles. If the slope becomes zero 
at the zero crossing of a singular thrust arc is encountered. This seems to be the case 

for Xf somewhat larger than 19 deg for the parameters chosen. Numerical difficulties are 

encountered which have to be dealt with in a special approach using higher order 
optimality conditions. This is beyond the scope of the present investigation. 
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The significant changes between the unconstrained and the heating constrained 
optimal synergetic maneuver are exhibited in the phase plots of Figure 9: While in the 
unconstrained case the heading angle varies almost linearly with the flight path angle over 
the central portion of the trajectory, in the constrained case the flight path angle is 
almost zero (Fig. 9a) there. The velocity is decreased and after thrust ignition increased 
again at almost constant altitude; this altitude is almost independent of the inclination 
change requested and is given grossly by the heating constraint for a lift coefficient 
correspond i ng to maximum L/D. In the unconstrained case the minimum altitude 
decreases with increasing plane change angle (Fig. 9b) and thrust is switched on at 
positive flight path angles. 

The optimal control time histories for the heating constrained skips are shown in 
Figures 8c through e for final heading angles of 7.5 to 19 deg. The plane changing thrust 
angle e is seen to be in the vicinity of tan' 1 (D/L) = 24.3 deg for L/D = 2.22 as has been 
estimated by the analytical approximate impulsive solution of References 19 and 20. 
There is an upward thrust component in the vertical plane ( a around 10 to 15 deg in 
Figure 8c) the average of which decreases with increasing plane change angle as has been 
indicated from the impulsive approximation [19,20], The aerodynamic bank 
angle p (Fig. 8d) is in the region of 60 to 85 deg for portions of the trajectory with 
relatively high dynamic pressure. The small initial values are due to the steeper entry 
angle chosen here as will be discussed later with reference to Figure 12b. The lift 
coefficient (Fig. 8e) varies around the value for (L/D) with characteristic jumps in 

ITlclX 

the slope as the heating constrained arc is entered or left. The high values of at the 

end occur when the vehicle has already left the denser atmosphere. The last corner in 
C^(t) is due to the end of the bum phase. In Figure 10 the mass ratio left after a plane 

change angle i is given for a complete maneuver to turn a 556 km (300 n.mi.) circular 
orbit. The exit velocity of the rocket engines on which the figure is based is 4.36 km/s. 
It is seen that the losses caused by the heating constraint are relatively small. The 
synergetic plane change maneuver is superior to an extra-atmospheric one for plane 
change angles i larger than 5 deg. For i = 20 deg the mass ratio in the new orbit is about 
15 percent higher after the heating constrained synergetic maneuver, which corresponds 
to a propellant reduction of about one-third for a shuttle-orbiter type vehicle compared 
to the extra-atmospheric, one-impulsive plane change. 


Sensitivity to Parameter Changes 

The most important aerodynamic vehicle parameters are the wing 
loading m Q /F at the beginning of the skip, the maximum L/D and the heating 

constraint. Another important parameter for the skip is the entry flight path angle 7 q . 
Their influence will be discussed for the reference skip trajectory given in Table 2. 

An increase in wing loading increases the ratio of inertial to aerodynamic forces 
for the same state and control of the vehicle. To decelerate the vertical velocity 
component to zero at the same altitude given by the heating constraint, the bank angle is 
reduced initially. At the same time, the switching function for the thrust is increased, 
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TABLE 2. REFERENCE SKIP TRAJECTORY PARAMETERS FOR 

THE SENSITIVITY ANALYSIS 


Wing Loading 

Drag Polar Coefficients 


m /F = 275 kg/m 2 


C Do = 04 ; k= 1 ; n= 1.86 


yielding (L/D) = 2.22 

: 6 v / ' max 


Heating Constraint Parameter AC^ = ^ ’ corres P on ding to 

T= 1093°C(2200°F) 


Jet Velocity 


Mass Flow Rate 


U = 4.36 km/s 
e 


m/m Q = 2 ' 1 O' 3 [s' 1 ] 



Initial 

Final 

Initial 

Adjoint Variable 

Final 

(Costate) 

Velocity, km/s 

7.95 a 

7.95 a 

-0.1656 

-0.1610 

Heading Angle, deg 

0 a 

15 a 

-0.516 

-0.636 

Flight Path Angle, deg 

-1.25 a 

2.0 a 

-0.289 

-0.266 

Lateral Range, deg 

0 a 

3.83 

0.374 

0 a 

Altitude, km 

80 a 

80 a 

-0.0002453 

-0.0003278 

Mass Ratio 

l a 

0.8/23 3 

-0.7939 

-l a 


a. Value is preset boundary conditions. 

leading to prolonged initial burn times for steep entry angles, thereby supporting the 
downward deceleration with a positive thrust angle a and reducing the wing loading for 
the aerodynamic flight phase. The total propulsion requirement increases slightly more 
than linearly with m Q , about one third of a percent drop in the final mass ratio for an 

increase in wing loading from 275 to 300 kg/m 2 . 
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The analytical approximate solutions for the unconstrained skip indicate that 
the AV requirement depends on the maximum lift to drag ratio by the 
factor exp(-D/L). A numerical test computation validates this dependency also for the 
heating constrained case. 


The effect of changing the severeness of the heating constraint is shown in Figure 
11. The adjustment parameter AC, h in equation (7) corresponds roughly to a limit 


temperature change. Letting AC 


LH 


= 0.03 is an approximation to increasing the limit 


temperature from 1039°C (2000° F) to about 1204°C (2200° F) for the given data. It 
shows how a limit temperature increase would affect the trajectory (Fig. 11a) and the 
controls (Figs, lib and 11c). The final mass fractions for the three-dimensional 
atmospheric/thrusting maneuver are 


4C lh 

0 

0.03 

0.06 

Unconstrained 

mJm 
P 0 

0.8123 

0.8145 

0.8162 

0.8200 


Due to the heating constraint with AC^^ = 0, the minimum altitude is lifted 

from 54 to about 61 km and the maneuver time increases from about 350 to 526 sec (50 
percent). The thrust switching function at t = 0 decreases for increasing AC^ , 

indicating that initial burn arcs become less favorable for hotter skips. At the same time, 
the initial bank angle increases (Fig. lib). From ix «= 90 deg in the unconstrained case, it 
is seen [ju = tan ' 1 (Xy/Xy C0S7)] that the entry flight path angle 7 = -1 .25 deg is 

Opt 0 

optimal for this parameter set since Xy Q « 0. The lift coefficient time histories again 
show the characteristic dips on the constrained arcs (Fig. 11c). 

The effect of changing the entry flight path angle is shown in Figure 12. The 

initial short bum phase appears only for entry angles steeper than approximately -1.28 

deg. As the entry angle becomes smaller, the ingoing leg of the maneuver changes while 

the outgoing leg remains essentially unchanged but is shifted in time. In the first half of 

the trajectory the bank angle increases to about 90 deg using almost all of the lift for 

plane changing. The maximum final mass fraction is obtained for 7 -0.78 deg and is 

c 

about 0.12 percent higher than for 7 =-1.25 deg. This result is obtained without 

considering viscous interaction effects. With these effects taken into account a steeper 
entry angle is advantageous [19, 20]. By considering the total heat transfer and the total 
maneuver time, an entry angle of -1 to -1.25 deg seems to be a good compromise. 
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CONCLUSIONS 


The numerical investigation of the synergetic plane change maneuver for a Space 
Shuttle orbiter-type vehicle under realistic constraints shows that, for the plane change of 
a 556-km (300-n.mi.) circular orbit over more than 5 deg, the synergetic maneuver 
requires less propulsion than an extra-atmospheric one. The kinetic heating constraint 
results in a markedly different trajectory shaping with much lower maximum dynamic 
pressure (one-third of the unconstrained case for Xf = 17deg) and normal acceleration 

(only 0.43 g^ compared to 1.6 unconstrained), but only a small decrease in the final 

mass fraction. 

The distribution of the burn arcs is very much dependent on the entry and exit 
flight path angles given as boundary values. For the nominal values chosen from total 
range and heat input considerations, an initial burn arc at entry may occur which 
decreases the flight path angle and reduces the wing loading for the aerodynamic 
maneuver. If the exit flight path angle is specified too low, the ascent burn phase may 
have two arcs, one when leaving the denser atmosphere (with an upward accelerating 
component) and a short second one at the right boundary. The change of the switching 
function for the thrust when approaching a critical constraint plane change angle, 
depending on boundary values and vehicle parameters, indicates that a singular thrust arc 
may be encountered. This question has to be examined in more detail. 

The results obtained for a nonrotating earth will be a good approximation for 
meridional entry in higher latitudes. For plane changes of 20 deg, a saving of 1 5 percent 
of the initial mass in propellant corresponding to a one-third propellant reduction may be 
achieved with the synergetic maneuver as compared to an extra-atmospheric one. 

It should be noted that the analysis in this document assumes classical point mass 
vehicle behavior. The effect of thrust and aerodynamic moments about the vehicle’s 
center of gravity should be examined in future studies to determine if it has any 
significant effect on the trajectory. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, March 1973. 
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Figure 1 . Synergetic plane change maneuver. 
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QUALITATIVE ALTITUDE CONSTRAINT FOR REENTRY 



BIVARIATE HEATING CONSTRAINT APPROXIMATIONS 


Figure 2. Kinetic heating constraint for Space Shuttle orbiter-type vehicle. 
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Figure 3. Coordinate system. 
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h/km 



FLIGHT PATH - HEADING ANGLE PHASE PLOT 



v/km/s 

ALTITUDE VELOCITY PHASE PLOT 

Figure 6. Multiple shooting optimization of the synergetic plane change. 

L/D = 2.22, xf = 35 deg, U e = 4.36 km/s, a Q ~ 1 g Q , v e = v f = 7.95 km/s, 

7 e = -1.25 deg, 7 f = 2 deg, and h e = h f = 95 km. 
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Figure 6. (Concluded). 



ALTITUDE TIME HISTORY FOR UNCONSTRAINED 
PLANE CHANGE WITH FINITE THRUST 

Figure 7. Unconstrained synergetic plane change with finite thrust (aj 0 
v e = v f = 7.95 km/s, y e = -1.25 deg, and m 0 /F = 275 kg/m 2 . 
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APPENDIX. COSTATE DIFFERENTIAL EQUATIONS 


From equations (18) and (19) follows for the differential equations for the 
Lagrangian multipliers (costate variables): 


; _ cos 7 
X v r~ 


cos x X x tan A - \ y 


('•*) 


- Xa sin x 


- X^ sin y 


+ -5^- e'^ h 

2 m 


( 2 VX v c D + C L w aj 


X^ = ~ V C ~ T ^X x sin x tan A + X A cos x 


) 



r 

— 

(x v p - x h v^ + y sin 7 

M 

1 - + X A sin X - X x cos x tan A 


tan e ^ h C L v si " * 


i _ V cos y COS x -V 

A A i . A-, 

A r cos 2 A X 


X h = '-p- ^X v sin 7 + cos 7 ^ + -p- cos 7 ^X y + X A sin X - X x cos x tan A^ 


V ( vX v C D + C L w a) 


X = _ e -0h v 

m 2m 2 


VX V C D + C L W 


On thrusting arcs, the following terms have to be added: 
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T V S ’ n 7 

X T = X^ + — 

T 1 7 mv cos 3 y Wr ^ 

X mT = ' ~^7* [( vX v) + (host) + X t 

Additive terms for the heating constrained arcs were given as equations (40), (4 1 ), and (42). 
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